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Abstract 

The Navier-Stokes-Coriolis system is a simple model for rotating fluids, which allows to 
study the influence of the Coriolis force on the dynamics of three-dimensional flows. In this 
paper, we consider the NSC system in an infinite three-dimensional layer delimited by two 
horizontal planes, with periodic boundary conditions in the vertical direction. If the angular 
velocity parameter is sufficiently large, depending on the initial data, we prove the existence 
of global, infinite-energy solutions with nonzero circulation number. We also show that these 
solutions converge toward two-dimensional Lamb-Osecn vortices as t — ► oo. 
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1 Introduction 



In recent years a lot of activity has been devoted to the mathematical study of geophysical flows, 
and in particular to various models of rotating fluids. Taking advantage of the stratification effect 
due to the Coriolis force, significant results have been obtained which are still out of reach for 
the usual Navier-Stokes system, such as global existence of solutions for large initial data [H [3] 
and stability of boundary layers for small viscosities [10\ I15j. We refer the interested reader to 
the recent monograph [4] which contains a general introduction to geophysical flows, an overview 
of the mathematical theory, and an extensive bibliography. 

In this article we study the so-called Navier-Stokes-Coriolis (NSC) system in a three-dimen- 
sional layer delimited by two infinite horizontal planes, assuming as usual that the rotation vector 
is constant and aligned with the vertical axis. This is a reasonably simple model for the motion 
of the ocean in a small geographic zone at mid-latitude, where the variation of the Coriolis force 
due to the curvature of Earth can be neglected. More realistic systems exist which take into 
account the variations of temperature and salinity inside the ocean, and include boundary effects 
modelling the influence of coasts, the topography of the bottom, or the action of the wind at 
the free surface, see (9j[I7]. Nevertheless, keeping only the Coriolis force is meaningful in a first 
approximation, because its effect is very important on the ocean's motion at a global scale due 
to the fast rotation of Earth compared to typical velocities in the ocean. 

Our main goal is to investigate the long-time behavior of the solutions to the NSC system 
for a fixed, but typically large, value of the rotation speed. As in [HE] we shall use the effect 
of the Coriolis force to prove global existence of solutions for large initial data, but the long- 
time asymptotics of those solutions turn out to be essentially two-dimensional and are therefore 
not affected by the rotation. Thus we shall recover as a leading term in our expansion the 
Lamb-Oseen vortex which plays a similar role for the usual Navier-Stokes system in the plane 
R 2 [8] or the three-dimensional layer R 2 x (0, 1) [19]. To avoid all problems related to boundary 
layers, we shall always assume that the fluid motion is periodic in the vertical direction. This 
hypothesis has no physical justification and is only a convenient mathematical way to disregard 
the influence of the boundaries. Although boundary conditions do play an important role in the 
problem we study and will have to be considered ultimately, in this paper we chose to focus on 
the motion of the fluid in the bulk. 

We thus consider the Navier-Stokes-Coriolis system in the three-dimensional layer B = R 2 x 
T 1 , where T 1 = R/Z is the one-dimensional torus. The points of B will be denoted by (x,z), 
where x = (x\,X2) £ R 2 is the horizontal variable and z E T 1 is the vertical coordinate. The 
system reads 

d t u + (u ■ V)u + Qe 3 A u = An - Vp , div u = , (1) 

where u = u(t, x, z) E R 3 is the velocity field of the fluid, and p = p(t, x, z) E R is the pressure 
field. Here and in what follows, it is understood that differential operators such as V or A 
act on all spatial variables (x,z), unless otherwise indicated. System ([T]) differs from the usual 
incompressible Navier-Stokes equations by the presence of the Coriolis term Qe^ A u, where 
fi E R is a parameter and = (0, 0, 1)' is the unit vector in the vertical direction. This term is 
due to the fact that our reference frame rotates with constant angular velocity 0/2 around the 
vertical axis. Note that ([I]) does not contain any centrifugal force, because this effect can be 
included in the pressure term —Vp. For simplicity, the kinematic viscosity of the fluid has been 
rescaled to 1, and the fluid density has been incorporated in the definition of the pressure p. 

As in the ordinary Navier-Stokes system, the role of the pressure in ([I]) is to enforce the 
incompressibility condition divu = 0. To eliminate the pressure, one can apply to both sides the 
Leray projector P, which is just the orthogonal projector in L 2 (B) 3 onto the space of divergence- 



2 



free vector fields. This operator has a rather simple expression in Fourier variables, which will 
be given in Appendix |Aj The projected equation then reads: 

d t u + P((« • V)u) + OT(e 3 A u) = Au , div-u = . (2) 

Another possibility is to consider the vorticity field u = curlu, which satisfies the following 
evolution equation: 

dtUJ + (u ■ V)oj — (oj ■ V)u — Qd z u = Auo . (3) 

Due to the incompressibility condition, the velocity field u can be reconstructed from the vorticity 
<jj using the Biot-Savart law, which in the domain B has also a simple expression, see Appendix lAl 
As is clear from (J3|), the vertical coordinate z plays a distinguished role in our problem 
because the rotation acts trivially on z-independent velocity fields. As a matter of fact, even if 
rotation is absent, the linear evolution dtu = Au leads to an exponential decay of the fluctuations 
of u in the vertical direction, due to the Poincare inequality. For these reasons, it is appropriate 
to decompose the velocity field as u(t, x, z) = u(t, x) + u(t, x, z), where 



u(t, x) = (Qu)(t, x) = I u(t,x, z) dz (4) 

is the average of u with respect to the vertical variable, and the remainder u = (1 — Q)u has 
zero vertical average. We shall say that u is a two-dimensional vector field in the sense that 
it depends only on the spatial variable x G M 2 , not on z, but one should keep in mind that 
u is not necessarily horizontal because its third component u^, is usually nonzero. A similar 
decomposition holds for the vorticity, and it is easy to verify that ui = curl it and Co = curlu. 
In particular, since d\U\ + 02^2 = and d\U2 — d^ui = ^3, the horizontal part of the two- 
dimensional velocity field u can be reconstructed from the third component of the vorticity tu 
via the two-dimensional Biot-Savart law, see Appendix^ This means that the averaged velocity 
field x) can be represented by two scalar quantities, namely us(t,x) and Cj^{t,x). 

We shall solve the Cauchy problem for equation ([2j) in the Banach space X defined by 

I={ae<(D) 3 divu = 0, u£H\B) 3 , usGtf^M 2 ), Cu 3 £ L 1 (R 2 ) n L 2 (R 2 )} , (5) 

equipped with the norm 

\\ U \\X = |H|#1(D) + 1 1^3 1 1 #1(^2) + ||(I> 3 || L 1( R 2) + 1 1 U 3 ||i2( R 2) . 

Observe that X <f_ ff 1 (D) 3 , because the two-dimensional horizontal velocity field Uh = (ui,U2) 
is not assumed to be square integrable. This slightly unusual choice is motivated by our desire 
to include infinite- energy solutions, which play a crucial role in the long-time asymptotics of the 
Navier-Stokes equations 0(8]. The most important example of such a solution is the Lamb-Oseen 
vortex, whose velocity and vorticity fields are given by the following expressions: 

u G (t,x) = -^=U G [-^=) , where U G (0 = - | 6 | , (6) 




wG (^) = —t G v~7m) ' where G(0 = ^ e "" r/4 1 1 • (7) 

As is easily verified, for any a£R and any Q G R, the vortex u(t, x, z) = au G (t, x) is an exact 
solution of the NSC system ©. In fact, one has F(u G ■ V)u G = and P(e 3 A u G ) = 0, so that 
u G solves the linear heat equation dtu = Ait. 

We are now in position to formulate our main result: 
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Theorem 1.1 For any initial data uq £ X , there exists Q,q > such that, for all £ R mt/i 
1^1 > ^ e iV5C system ([2]) /ias a unique global (mild) solution u £ C°([0, oo),X) satisfying 
u(0) = uq. Moreover \\u(t, •) — au G (t, -)\\x — > as t — > oo, where 

a = (cuTluo)sdxdz . (8) 

This theorem contains in fact two different statements. The first one is the existence of global 
strong solutions to the NSC system (|2|) for arbitrarily large initial data in X, provided that the 
rotation speed |0| is sufficiently large (depending on the data). To prove this, we closely follow 
the existence results that have been established for rotating fluids in the whole space M 3 , see [U 
Chapter 5]. In particular, if the three-dimensional part u of the solution is not small at initial 
time, we assume that the rotation speed |0| is large enough so that u is rapidly damped by the 
dispersive effect of the linearized equation 

d t u + OP(e 3 A u) = An , divu = . (9) 

For the reader's convenience, we briefly recall in Section 12.21 and Appendix [B] the Strichartz 
estimates satisfied by the solutions of ([9]) with compact support in Fourier space. Except for the 
choice of the spatial domain, the main difference of our approach with respect to [4] is that we 
do not assume that the whole velocity field u belongs to L 2 (B) 3 . As a consequence, we cannot 
use the energy inequality which plays an important role in the classical approach. To guarantee 
that the two-dimensional Navier-Stokes system has uniformly bounded solutions, the hypothesis 
u~h = (1*1,1*2)* £ L 2 (R 2 ) 2 is replaced by 0)3 £ L 1 (IR 2 ), a condition which allows for solutions with 
nonzero total circulation such as the Oseen vortex (JSJ) , (JZJ) • 

The second part of Theorem 11.11 which concerns the long-time behavior of the solutions, is 
more in the spirit of the previous works [HJ [19] . When stated more explicitly, our result shows 
that the solution u(t, x, z) satisfies 



lifl(B) + IP3W||hi(IR2) + 11^3(^)11^2^2) — — * 

and 

i+t Vyr+I 



u 3 (t)--^-g[ . ) >0, (10) 



where g(£) = (47r) _1 e / 4 . In particular, if the total circulation a is nonzero, we see that 
0)3 (t) does not converge to zero in the (scale invariant) space L 1 (M 2 ), but to the Oseen vortex 
with circulation a, which is thus the leading term in the asymptotic expansion of the solution 
as t — » 00. This is in contrast with the case of finite-energy solutions, which always converge to 
zero in the energy norm. 



We conclude this introduction with a few additional remarks on the scope of Theorem 11.11 

1) As is well-known, it is possible to prove the existence of solutions to the NSC system (J2|) 
under weaker assumptions on the initial data. For instance, it is sufficient to suppose that 
u(0) £ F 1 / 2 (D) 3 , u 3 (0) £ L 2 (M 2 ), and a> 3 (0) £ L 1 (M 2 ), in which case the solution u(t) will 
belong to X for any positive time. Since we are mainly interested in the long-time behavior of 
the solutions, we disregard these technical details and prefer working directly in the (noncritical) 
space X. 

2) Theorem 1 1 . 1 1 does not give any information on the convergence rate towards Oseen's vortex. 
The proof shows that ||V-u 3 (i)|| L 2( K 2) + ||w3(*)IIl 2 (r 2 ) = 0(t~ 1 / 2 ) and ||«(i)||#i(B) = 0(e~ vt ) for 
all v < 4-7T 2 as t — ► 00, but without additional assumptions on the data it is impossible to specify 
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the decay rate of ||^3(t)||L 2 (R 2 ) or the convergence rate in (fTU|) . However, algebraic convergence 
rates can be obtained if we assume that the initial data uq(x) decay sufficiently fast as \x\ — ► oo, 
see [HI US]. 

3) In the proof of Theorem 11.11 we need a large rotation speed only to prove the existence 
of a global solution, in the case where uq = (1 — Q)uq is not small. Once existence has been 
established, the convergence to Oseen's vortex holds for any value of O and does not rely on the 
Coriolis force at all. Since our domain B has finite extension in the vertical direction, we can 
use the Poincare inequality to show that u{t) converges exponentially to zero as t — > oo, but 
this point is not crucial: Our proof can be adapted to cover the case of the whole space M 3 , if 
we assume as in [3] that u = u + u with u G i7 1 (M 3 ), or even u G i7 1/,2 (IR 3 ). In this situation 
the decay of u(t) will not be exponential. 

4) As is explained in [19], we can prove the analog of Theorem Oin the layer M? x (0, 1) with 
different bounday conditions, for instance stress-free conditions. The case of no-slip (Dirichlet) 
boundary conditions is very different, because the solutions will converge exponentially to zero 
as t — > oo, and the Oseen vortices can only appear as long-time transients. 

5) A careful examination of the proof shows that the angular velocity Oo i n Theorem 11.11 can 
be chosen in the following way: 

O = max^||V£t ||i2 - K , o) , with K = Ce c||uo11 * , 

where uq = (1 — Q)uq and C > is a universal constant. In particular, one can take Oo = if 
uq is sufficiently small, depending on uq. Of course, there is no reason to believe that this result 
is sharp. 

The rest of this paper is organized as follows. In Section [2] we prove the existence part of 
Theorem 11.11 using energy estimates for the full system ([2]) and dispersive (Strichartz) estimates 
for the Rossby equation @ . Section [3] is devoted to the convergence proof, which relies on 
a compactness argument and a transformation into self-similar variables. In Appendix |A] we 
collect a few basic results concerning the Biot-Savart law in the domain D, and in Appendix [Bl 
we give a proof of the dispersive estimates for equation ([9]) which are used in the global existence 
proof. 

Acknowledgements. The authors are indebted to Isabelle Gallagher for helpful discussions 
on several aspects of this work. 

2 The Cauchy problem for the Navier-Stokes-Coriolis equation 

In this section we prove that the Navier-Stokes-Coriolis system ([2]) is globally well-posed in 
the function space X defined by ©, provided that the rotation speed O is sufficiently large 
depending on the initial data. The precise statement is: 

Theorem 2.1 For any initial data uq G X, there exists Oo > such that, for all O G R 
with |0| > Oo, the NSC system ([2|) has a unique global solution u G C°([0, oo),A) satisfying 
u(0) = uq. Moreover, there exists C > (depending on uq) such that \\u(t)\\x < C for allt > 0. 

As is clear from the proof, one can take Oo = in Theorem 12.11 (hence also in Theorem II. lj) 
if the three-dimensional part no = (1 — Q)uq of the initial velocity field is sufficiently small in X, 
see Remark 12.101 below. For large data, however, nobody knows how to prove global existence 
without assuming that the rotation speed O is large too. 
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2.1 Reformulation of the problem 



If u(t, x, z) is any solution of the NSC system ([2]), we decompose 

u(t,x,z) = u(t,x)+u(t,x,z) , (11) 

where u = Qu, u = (1 — Q)u, and Q is the vertical average operator defined in (jlj). Our first task 
is to derive evolution equations for u and u. Integrating ([2]) over the vertical variable zGT 1 , 
and using the fact that P and Q commute with each other (see Appendix |Aj), we obtain 

d t u + ¥[(u- V)m + Q(m- V)fi] = Am , divn = . (12) 

This is a two-dimensional Navier-Stokes equation for the three-component velocity field u(t,x), 
with a quadratic "source term" depending on u. Remark that the Coriolis force disappeared 
from (|12p . because curl(e3 Am) = —d z u = 0, so that P(e3 Am) = 0. On the other hand, subtracting 
(fTD from ([2]), we find 

d t u + ¥[(u- V)m + (m- V)m+(1-Q)(m- V)m] + OP(e 3 Am) = Am , divM = . (13) 

Thus m(£, x, z) satisfies a three-dimensional Navier-Stokes-Coriolis system, which is linearly cou- 
pled to (|12|) through the transport term P(m • V)m and the stretching term P(m • V)m. 

As is explained in the introdution, the averaged velocity field u(t, x) can be represented by 
two scalar quantities, namely its vertical component u~ 3 (t,x) and the third component Q 3 (t,x) 
of the averaged vorticity field. Taking the third component of (fT2"|) and using the fact that 
(Pm)3 = u 3 (see Appendix |Aj), we obtain the following evolution equation: 

d t u 3 + (u h ■ V)m 3 + Ni = Au 3 , xeM 2 , t > , (14) 

where Uh = (ui,U2Y and Aq = Q(m • V)m3. Similarly, if we take the third component of ([3]) and 
integrate the resulting equation over the vertical variable z, we find 

d t 3 + (u h ■ V)u> 3 + N 2 = Auj 3 , i£K 2 , t > , (15) 

where = Q((m • V)lD3 — (Co ■ V)u 3 ). Here we have used the fact that (u) ■ V)u 3 = 0, see (jT9|) 
below. 

By construction, the original NSC equation ([2]) is completely equivalent to the coupled 
system (fl~3|) . (HU), (fT5|) . To prove local existence of solutions, we consider the integral equations 
associated to these three PDE's (via Duhamel's formula), and we apply a standard fixed point 
argument in the function space C°([0, T], X). The result is: 

Proposition 2.2 For any r > 0, there exists T = T{r) > such that, for any H 6 1 and all 
initial data uq S X with ||mo||x < r , the Navier-Stokes-Coriolis system §2§ has a unique local 
solution u £ C°([0, T],X) satisfying m(0) = Mo- 

The proof of this statement uses classical arguments, which can be found in [6], [11], |12j . 
and will therefore be omitted here. The fact that the local existence time T depends on mo only 
through (an upper bound of) the norm ||mo||x is not surprising, because we work in a function 
space X which is not critical with respect to the scaling of the Navier-Stokes equation. However, 
it is worth noticing that T is independent of the rotation speed Q. This is because the rotation 
does not act at all on the two-dimensional part (|14|) . (|15|) of our system, whereas in (|13|) it 
appears only in the term W(e 3 Am), which is skew- symmetric in the space ^(B) 3 and therefore 
does not affect the estimates. 
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To prove global existence and conclude the proof of Theorem l2,ll it remains to show that any 
solution u G C°([0,T], X) of (|2j) is bounded for all t G [0, T] by a constant depending only on 
the initial data no = u(0). As is well-known, this is relatively easy to do if the three-dimensional 
part no of the initial data is small in ^T 1 (D), see [6], [13) . In the general case, we shall use the 
dispersive properties of the Rossby equation ([9]) to prove that the solution u(t,x, z) of (fl~3j) is 
rapidly damped for positive times if the rotation speed is sufficiently large. 



2.2 Dispersive properties 

Since our spatial domain B = R 2 x T 1 is bounded in the vertical direction, the Poincare inequality 
implies that the solutions of the linear equation ([9]) decay exponentially to zero as t — > do. More 
precisely, for any s > and all divergence-free initial data uq G (1 — Q)H S (B)) 5 , the solution 
u(t, x, z) of © satisfies 

\\u(t)\\ H s m < \\uq\\h'(D) e" 4 " 2 * , t > . (16) 

This estimate is straightforward to establish by computing the time-derivative of ||n(t)||^ s and 
using the Poincare inequality ||Vn||^ s > 47r 2 ||u|| 2 ^ s together with the fact that the Coriolis 
operator u t—* P(e3 A u) is skew-symmetric in iP(B) 3 for divergence-free vector fields. Note 
in particular that (I16p is independent ot the rotation speed CI. However, as is shown e.g. in 
[3], additional information can be obtained for large \Cl\ if we exploit the dispersive effect of 
the skew-symmetric term n¥(e^ A u). The corresponding Strichartz-type estimates are most 
conveniently derived if we restrict ourselves to solutions with compact support in Fourier space. 

Throughout this paper, we use the following conventions for Fourier transforms. If / G L 2 (B) 
or L 2 (B) 3 , we set 

f(x, z) = ±- ! Y, fn{k) e <(*-*+*"«) dk , x G R 2 , z G T 1 , (17) 

where 

f n (k) = J- / / f(x,z)e~ i( - k - x+2nnz ^ dzdx , k£R 2 , nGZ. (18) 
2vr J R 2 Jji 

With these notations, the norm of / in the Sobolev space H S (B) can be defined as 

\\f\\ H s = (f ^(1 + |A;| 2 + 47rVr|/ n (A;)| 2 dA : ) 1/2 , (19) 

where \k\ 2 = k\ + k 2 . Given any R > 0, we denote by Br the ball 

Br = { (k, n) G R 2 x Z y/\k\ 2 + 47r 2 n 2 < r} . (20) 
Following closely the approach of [4, Chap. 5], we obtain our main dispersion estimate: 



Proposition 2.3 For any R > 0, there exists Cr, > such that, for all uq G (1 — Q)L 2 { 
with divno = and supp (uo) n (k) C Br, the solution no/© with initial data no satisfies 



^•W L 1 (WL+ ,L°° (B)) - Cr\CI\ 4 \\u \\ L 2, B ) . (21) 



For completeness, the proof of this proposition will be given in Appendix [Bl Estimate (|2ip 
clearly demonstrates the dispersive effect of the Coriolis term in @: If the initial data no are 
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compactly supported in Fourier space, the L°° norm of the solution u(t, •) will be very small (for 
most values of time) if the rotation speed |0| is large enough. This is in sharp contrast with 
what happens for Sobolev norms, for which the best we can have is estimate (|16p . As a side 
remark, if we consider initial data no whose Fourier transform is supported outside the ball Br, 
then we clearly have 

IK*)||ff« < WMIh* e~ R 1 for all t > 0. 
Combining Proposition 12.31 with estimate ()16p . we deduce the following useful corollary: 

Corollary 2.4 Under the assumptions of Proposition \2.'J\ the solution u of © satisfies, for 
any p G [1, +oo] and any q G [2, +oo] such that — h ~ < 1, 

i_ 

II«IIlp(R+,L9(d)) < C R (Q) 4 " ||n || L 2 (]D) ) , (22) 

where (ty = (l + l^ 2 ) 1 / 2 . 

Proof. Fix s > 3/2. Using Sobolev's embedding and our assumptions on uq, we obtain from 

< C\\u(t)\\ H s < C\\u \\ H s e~^ H < C R \\u \\ L 2 e~^ H , t>0, (23) 

where Cr denotes a generic positive constant depending only on R. In particular, we have the 
estimate ||u||li(k < Cr||^o||l 2 fo r all £ 1, so that OTj) holds with |f2| replaced by (O). 
This gives ([22]) for (p, q) = (l,oo), and since the case (p, q) = (oo,oo) is immediate from ([23]) . 
we see that (|22l) holds for all p G [l,oo] if q = oo. Finally, as ||n(t)||i2 < ||no||L 2 f° r au t, > 0, 
the general case follows by a simple interpolation argument. ■ 

To exploit the dispersive properties of the linear equation ([9]) in the analysis of the nonlinear 
problem (|13p . we use the following decomposition, which is again borrowed from [1]. Let x £ 
C^°(1R) be a cut-off function satisfying < x( x ) < 1 fo r an x G M, x( x ) = 1 fo r \ x \ < 1/2 and 
x(x) = for |x| > 1. Given any R > 0, we define the Fourier multiplyer Vr = x(|V|/i?) by the 
formula 

(W)=x(^^)/nW, ^GM 2 , (24) 
If u(t, x, z) is a solution of (|13p with initial data «o(x, z), we decompose 

n(t,x,z) = X(t,x,z) +r(t,x,z) , (25) 

where A(t, x, z) satisfies the linear Rossby equation 

d t X + nP(e 3 A A) = AA , divA = , (26) 

with initial data Ao = Vru$. By construction, the remainder r(t,x,z) is a solution of the 
nonlinear equation 

d t r + QF(e 3 Ar) + N 3 = Ar , divr = , (27) 

with initial data r = (1 - Vr)u , where N 3 = ¥[(u ■ V)n + (u ■ V)u + {l-Q)(u- V)u]. 

In the rest of this section, we consider equations (|26p . (I27p instead of (|13fl . so that our final 
evolution system consists of (fl"4"j) . (fT5|) . (p6|) . (p7|) . Given uq = uq + tto G X, we will choose the 
parameter > large enough so that the initial data ro = (1 — Vr)uo for equation (f27|) are 
small in ff 1 (D). Then the rotation speed will be taken large enough so that we can exploit 
the dispersive estimates for X(t,x,z) given by Corollary 12.41 
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2.3 Energy estimates 

We now derive the energy estimates which will be used to control the solutions of the nonlinear 
equations CH]), ([IS]), flZT 



Proposition 2.5 There exist positive constants Cq,C\ such that, if u 6 C°([0,T], X) is a solu- 
tion of © for some H 6 K, and if u is decomposed as in (fTT|) . ([25]) /or some R > 0, then the 
corresponding solutions of (THj) . (fT5|) . (f27|) satisfy, for any t £ (0,T]: 



<l "n 3 (t)||i 2(IR2) < -||Vu 3 (t)||| 2(K2) + ||n(t)||l 4(B) , (28) 



di 1 
d? 

' IMt)Hl2 (R2) < -||Vo)3(t)||| 2(R2) + 8|p(t)||Vn(t)|||| 2(D) , (30) 



(1 , |V^ 3 (i)|li 2(M2) < -||A^3(t)|li2 (M2) + Co(||V^3(Olli 2 ||^3(i)||| 2 + |||^(t)||Vn(^|||2 2 ) , (29) 
'1 — ,2 



dt 1 



d 

dt 1 



I^3(*)||li(r2) < ||W3(0)|Ui(Ra) + 2 f ||tt(s)|| L 2 (D) ||Au(s)|| L 2 (D) ds , (31) 

J 

Vr(i)|| 2 2(D) < -||Ar(i)|| 2 2(D) +C 1 ||Vr(i)||| 2 ||V«(t)||i 2 ||An(t)||^ 2 (32) 

+ cx(ii«(i)ii| 4 ||vA(t)ni 4 + \\vm\hm)\\u + \\\m\ \^m\\\h) ■ 



Remark 2.6 Here and in what follows, if f is a vector valued or matrix valued function, we 
denote by \ f\ the scalar function obtained by taking the Euclidean norm of the entries of f . 
Given any p G [1, oo], we define as |||/|||lp. With these conventions, if to = curl-u, we 

have for instance \u\ < \/2\Vu\ and \\uj\\l 2 = \\^ u \\l 2 - 

Proof. To prove (I28p . we multiply both sides of (|14H by us and integrate over R 2 . The transport 
term {u^ ■ V)«3 gives no contribution, because u is divergence-free, and the diffusion term Ait3 
produces the negative contribution — ||Vn3||^ 2 after integrating by parts. Since 

— / U3N1 dx = — J us(u ■ V)tt3 dx dz = / n 3 (tt • VU3) dx dz < -\\u\\ 4 L 4 + -\\Vu3\\ 2 L2 , 
Jm 2 Jd Jd 2 2 

we obtain the desired estimate. In a similar way, to prove (f2"9"j) . we multiply (TH|) by — Au^ 
and integrate over R 2 . The transport term gives here a nontrivial contribution which, after 
integrating by parts, can be bounded as follows: 

/ (Au 3 )(u h ■ V)tt 3 dx < / |Vu 3 ||Vu/j||V?23|dx < ||Vn 3 ||| 4 ||Vu fe || L 2 
Jm 2 Jk 2 

< C||Att 3 || i2 ||Vit3||£2||a;3|| i2 < -||Aw 3 || 2 2 + C||Vw 3 ||| 2 1|^3|||2 • 

Here, to get from the first to the second line, we have used an interpolation inequality and the 
fact that Uh is obtained from Q3 via the Biot-Savart law (I8ip . see Appendix |A] Since we also 
have 

/ A-U3 N\ dx < / |An3||u||Vn| dxdz < -HA-usH^ + |||n||Vn|||^ 2 , 
Jm 2 Jd 4 

we obtain again the desired inequality. 

On the other hand, multiplying (TT5|) by 0)3 and integrating over R 2 , we easily obtain (fHUj) . 

because 

- / u> 3 N 2 dx = / (a) 3 (n • V)(D 3 - u 3 (d) • V)(D 3 ) dxdz < -||Vu)3||| 2 + 2|||u||u)||| 2 2 , 
Jr 2 Jd^ / 2 
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and \Cj\ 2 < 2|V-u| 2 . The prove (|3ip we observe that, since the vector field Uh is divergence-free, 
any solution of (fl~5j) in L 1 (M 2 ) satisfies 

\\u 3 (t)\\iA < IN(0)|| x i+ / ||^ 2 (s)|| r idfl, t>0. 

J o 

This bound can be established using the properties of the fundamental solution of the linear 
convection-diffusion equation dtf + (uh ■ V)/ = A/, which will be recalled in Section [3T2l below. 
Since 

H-^Hl 1 — II^IU 2 II^^IIl 2 + II^IIl 2 II^^IIl 2 — II^IU 2 l|A£t||x,2 + II^^IIl 2 — 2||w||i,2 1| Au 1 1 £,2 , 

we obtain (j3"Tj) . 

Finally, to prove (|32|) . we multiply (I27p with — Ar and integrate over B. As was already 
explained, the Coriolis term iW(e^ A r) gives no contribution, because it is skew-symmetric in 
any Sobolev space. So we just have to bound the contributions of the nonlinear term N3, which 
are threefold. Since u = A + r, the transport part ¥(u ■ V)u in N3 produces two terms, which 
can be estimated as follows: 



/ Ar • (u • V)Adxdz < — 1| Ar|| 2 2 + C\\u\\ 2 L4 \\VX\\ 2 L 4 , 
Jo 10 

/ Ar • (u ■ V)r dx dz < / |Vr| |Vu| |Vr| dx dz < ||Vr|| 2 8 ||V 
Jo Jo L3 



1 



< C||Vr||^||Ar||| 2 ||Vu||2 2 ||A«||£ 2 < C||Vr||£ 2 ||Ar||£ 2 ||Vu||£ 2 ||Au||| 2 
1 



< 



10' 



Ar||2 2 + CUVrllfallV-ulllallAu 



Il 2 I 



II 2 ■ 



Here we have used interpolation inequalities, Sobolev embeddings, and the Poincare inequality 
||Vr||£2 < C||Ar||j^a. The two terms produced by the stretching part ¥{u ■ V)u in N3 can be 
estimated in a similar way: 



/ 



Ar • (A • V)udxdz < — || Ar||| 2 + CHAHf^ ||Vu|| 2 2 



Ar • (r • V)u dx dz 



< || Ar 11^2 ||r||£4 ||Vw||£,4 < C|| Ar|| i2 1| Vr|| i2 ||Vu||? 2 II Au||^ 2 



< -||Ar||£ 2 +C||Vr||£ 2 ||Vu||£ 2 ||AzZ||£ 2 



Finally, the contribution of the quadratic term P(l — Q)(u ■ V)tt in N3 satisfies 



1 



/ Ar • (u • V)£tdxdz 
Jo 

Collecting all these estimates, we obtain (|32p . This concludes the proof. 



< — ||Ar[|£a +C[||fi| |Vw| 



2 

L- 



2.4 Global existence 

In this section, we combine the dispersive properties of Section 12.21 and the energy estimates of 
Section 12.31 to complete the proof of Theorem 12.11 We start with a preliminary result, which 
summarizes in a convenient way four of the five inequalities established in Proposition 12.51 
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Lemma 2.7 There exist positive constants C 2 , C3, and C4 such that the following holds. Let 
u G C°([0,T],X) be a solution of © for some Q G M., which is decomposed as in (|llj) . (I25D for 
some R > 0. Assume moreover that there exist K > 1 and e G (0, 1] smc/i i/iat t/ie corresponding 
solutions of (fT4"|) . (fT5|) . (12T1) safe/y 

||V« 3 (t)||r»CR a ) ^ ^ 2 . l|w3(*)||za(R3) < # , l|Vr(i)|| L2(D) < e , (33) 
/or i G [0, T] . // me define 

*(*) = ll«3(0lli2(K») + ll^3(i)lli 2( R2 ) + <5||Vn 3 ||i 2(K2) + ||Vr(t)||2 2(D) , (34) 
/or some (5 G (0, 1], i/ien 

^*(t) < -(||Vn3(t)||| 2 + ||VcD3(t)||| 2 +5||An 3 (t)||l 2 + ||Ar(t)||| 2 ) 

+ C ^ 2 ||Vn 3 (t)||i 2 + C 2 e 2 i^ 4 ||A U (t)||| 2 + C 3 e 2 1| Ar(t) ||| 2 (35) 
+ (5- l <S>{t)+K\\^{t)\\ L ,)G(t)+F{t)+e 2 G(t) , 

for all t G (0,T], w/tere 

F(t) = C 4 (||A(t)|| 4 L4 + ||A(t)||£ BO ||VA(t)||£ a ) , (36) 
G(t) = C 4 (||VA(t)||£. + ||A(t)||i- + ||VA(t)||| 4 ) . 

Proof. If <I> is defined by (|34|) . it follows immediately from Proposition 12.51 that 

^*(t) < -(||V%(t)||| 2 + ||Va)3(i)||i 2 +^||Au3(i)||| 2 + ||Ar(t)||| 2 ) 

+ ||n(t)||l 4 + C|||n(t)||Vn(t)|||i 2 +Co<5||Vn3(t)||l 2 ||^ 3 (Olli 2 (37) 
+ C(||Vr(t)||i a ||Vfl(t)||i a ||Au(t)||» a + ||u(t)|&||VA(t)||i4 + ||V<Z(i)||£ 2 .||A(*)ll£-) . 
Using interpolation inequalities, Sobolev embeddings, and the a priori bounds (|33p. we first get 

||n(t)|| 4 L4 < C(||r(t)|| 4 L4 + ||A(t)|| 4 L4 ) < C(||r(t)||| 6 ||r(t)|| L2 + ||A(t)||l 4 ) 

< C(||Vr(t)||| 2 ||r(t)|| L2 + ||A(t)|| 4 4 ) < Ce 2 ||Vr(t)|| L2 ||r(t)|| L2 + F 1 (t) , 

where F\(t) = C||A(i)||* 4 . Proceeding in the same way, we also obtain 

\\\n\\Vu\\\h < G(|||r||Vr|||| 2 + || \r\ | VA| ||| 2 + || |A| |Vr| ||| 2 + || |A| |VA| ||| 2 ) 

< C(||r||| 6 ||Vr||| 3 + llrl^llVAHloc + || Vr||| 2 [| A]||oo + ||A||ioc[|VA||^) , 
so that |||tt(t)||Vw(t)|||| a < Ce 2 ||Vr(t)|| L2 ||Ar(i)|| L2 + e 2 G 1 (t) + F 2 (t), where 

G 1 (t) = C(||VA(t)||i» + ||A(i)||L«) , F 2 (t) = C\\\(t)\\U\V\(t)\\l> ■ 

It remains to estimate the last four terms in the right-hand side of (|37|) . The first two in this 
group are independent of A, and are simply bounded using assumption (j33H and the fact that 
||V«||? 2 = ||V^3|lx,2 + ||^3|| 2 2 - On the other hand, in view of Proposition I A. H we have 

||u|| 2 4 < C(||u 3 ||| 4 + ||^3|| 2 4) < C(||""3 Hi, 2 II V r €t 3 ||z,2 + ||lD 3 || l i||o) 3 || l2 ) , 

Li 3 

hence 

||n(i)||l 4 ||VA(t)||| 4 < (5-^) + K\\u, 3 (t)\\ L1 )G 2 (t) , 

where G 2 (t) = C||VA(i)|| 2 4 . Similarly, we find ||W(t)||| 2 ||A(t)|||oo < S^^Gity). Thus, 
using the Poincare inequality ||r|| i2 < ||Vr||^ 2 < ||Ar|| i2 , we see that (I35p holds with F(t) = 
Fi(t) + F 2 (t) and G(t) = G x {t) + G 2 (t). ■ 
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Remark 2.8 In view of Corollary \2.4\ there exists a constant Cr > (depending only on R) 
such that 

/ F(t)dt < C R {n)~^\\u \\ 4 L2 , and / G(t)dt < C R {n)~^ \\u \\ 2 L2 . (38) 
Jo Jo 

Remark 2.9 Without loss of generality, we shall assume henceforth that the constants which 
appear in Proposition \2.5\ and Lemma \2. 7| satisfy Ci > 1, i = 0, . . . , 4. 



Proof of theorem 12.11 Given uq G X, we define u(0) = Quq, uq = (1 — Q)uq, and uj(0) = 
curl(Quo), where Q is the vertical average operator (jj]). We first choose K > 1 such that 

K 2 

II "3 (0) Ulrica) + ||W3(0)|| Z 2 (R2) + ||w 3 (0)|| L i(K2) + 2||«o||^i( D ) < J^T > ( 39 ) 



where Co > 1 is as in Proposition 12.51 Next, we take e G (0, 1] sufficiently small so that 

£2 * mm { 2^ ' 2C^I ' Wh6re 6 = 2G^ G ( °' 1] ' (40) 



,2"^'2C^}' Wh6re 5= 2C7^ 

and C2 > 1, C3 > 1 are as in Lemma f2,7l Once this is done, we set Ao = VrUq and ro = 
(1 — Vr)uo, where Vr is the Fourier localization operator defined by (j24"|) . We assume that the 
parameter R > is sufficiently large so that 



4e 2ClX8 ||Vr ||i 2 < e 2 , (41) 

and we denote by X(t,x,z) the solution of (I26D with initial data Ao- Finally, using Remark 12.81 
we choose f2o > sufficiently large so that, if > £lo, 

roo roo jf2 

/ G(t)dt < 51og(2) , / (F(t) + e 2 G(t))dt < -— , (42) 
4 e 2Cl ^ 8 J (F(t) + (K 4 + e 2 )G(t)) dt < e 2 . (43) 



and 



Remark 2.10 If uq is small enough so that 4e c ' 1 ^' 8 ||Vtto|li2 < £ 2 , then we can take formally 
R = 0, so that ro = u$ and Ao = 0. In that case, one has F(t) = G{t) = 0, and (I42|) . (|43l) are 
0/ course satisfied for any DgK. 

By Proposition 12.21 equation ([2]) has a unique maximal solution u G C°([0, T*), X) with 
initial data no, where T* G (0, +00] denotes the maximal existence time. If we decompose 
u(t) = u(t) + X(t)+r(t) as in (HU), ([25]), then <S 3 (i), u) 3 (i), r(t) are solutions of CEU), CE]), (l27l) . 
respectively, and we know from (|39p and (|41|) that 

K if e 

l|Vn 3 (0)|| L 2(R2) < -J , 11^(0)1^2^2) < — , ||Vr || L 2 (ID)) < - . 

Thus, by continuity, the bounds (|33H will be satisfied at least for t > sufficiently small. Let 



T = sup< T G [0,T*) 



The bounds ([Ml) hold for all t G [0,T]| G (0,T*] . (44) 



We shall prove that T = T*. This implies of course that T = T* = +00, and that the solution 
u(t) of ([2]) stays bounded in X for all t > 0, as is claimed in Theorem 12 .11 
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Assume on the contrary that < T < T*, and let ^(t) = $(i) + \\u)3(t)\\ L i, where $ is defined 
in ([Ml). Using (J55D and we find 

Ht) + i 7^ £ (H^^3(^)|li 2 ||V^ 3 (^)|li 2 ^||^^3(^)|li 2 + ||^( S )||i 2 ) (45) 

< $(0) + ^ 1 I ^{s)G{s)ds + I (F(s) + e 2 G(s))ds , te[0,T]. 
Jo Jo 

On the other hand, since 

2\\u{t)\\v\\*m\» < 2^l|An(t)||| 2 < ^(||Ar(t)||| 2 + ||AA(t)||| 2 ) , 
it follows from (1311) that 



||0 3 (t)|Ui < \\*3(P)\\i> + \J l|Ar( S )||| 2 d S + ||Vn |li 2 , *e [0,T] . (46) 

Here we have used the fact that 2 J °° ||AA(i)||| 2 dt = ||VA || 2 2 < ||Vtt ||^ 2 by ((26]). Summing 
up gSD and we obtain for t E [0, T): 

*(*) + i^*(" v ^c S )[ii 2 -h iivcD 3 ( s )iii 2 + ^[i^ 3 (^)[ii 2 -h |[|^( s )iii 2 ) a s (47) 

< ^(Oj + IIVfiollia / V(s)G(s)ds+ f (F(s) + e 2 G(s))ds . 



This integral inequality for ^(t) can be integrated using Gronwall's lemma. In view of ((39 
[I]) and (I42p . we easily obtain 



*(*) + i ^(HV^WII^ -+- [|Va3ft(»)]|i» H- tfHAfiaWHisi + §II^WIIi-) (48) 

< 2(*(0) + ||Vfi ||i a + j\F( S ) + e 2 G(s))ds) < ^ , 

for all t G [0, T]. In a similar way, using ((32j) . (f33l) and proceeding as in the proof of Lemma [2~7l 
we find 

l|Vr(t)||l 2 + ±£\\Ar(s)f L2 ds < ||Vr ||| 2 + 2C l K i J* ||Vr( S )||| 2 ||Au( S )||| 2 ds (49) 

+ (T 1 / tf(s)G(s)ds + / + e 2 G(s))ds , te[0,T]. 

io Jo 

From we know that Jj ||A«(s)||| 2 ds < 2K 2 /(4C 5) = ET 4 . Thus we can apply Gronwall's 
lemma to (I49D and, using in addition (I4ip and (I43p . we obtain 

l|Vr(i)|| 2 2 + i^||Ar( S )||| 2 d S (50) 
< e 2^^ 8 (||Vr ||| 2 +^ 4 ^G( S )d S + ^(F( S ) + £ 2 G( S ))d S ) < y , 

for all t € [0,T]. 

Now, it follows immediately from (1480 . (I50p that 

HViiaWHi. < ^ , < > llVr^lli, < y , 

for all i E [0,^], which obviously contradicts the definition (|44j) of T. Thus T = T* = +oo, and 
estimates ((53" ]) . (|15]> . ([50]) hold for all t > 0. This concludes the proof of Theorem EU ■ 
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3 Convergence to Oseen Vortices 



To complete the proof of Theorem 11.1] it remains to show that the global solution u(t, x, z) 
of the Navier-Stokes-Coriolis system ([2]) constructed in Section [2] converges to Oseen's vortex 
as t — > 00. To do that, we decompose u(t, x,z) = u(t,x) + u(t, x,z) as in (fTT|) . and we first 
show that the three-dimensional part u(t) converges exponentially to zero in H 1 ^) 3 , due to 
Poincare's inequality. We next turn our attention to the two-dimensional part u, and prove that 
the third component u 3 (t) decays to zero in ff 1 (M 2 ). Finally, the most delicate point is to show 
that Co 3 (t) converges to Oseen's vortex in L 1 (R 2 ) as t — > 00. Here the main ingredients are a 
transformation into self-similar variables, a compactness estimate for the rescaled solution, and 
a characterization of the complete trajectories of the two-dimensional Navier-Stokes equation 
which was obtained in [8]. 



3.1 Exponential decay of u 

We recall from ([23]) that u(t,x,z) = r(t,x,z) + X(t,x,z), where A satisfies the linear equation 
(|26D and r is a solution of ((27]). We already know that ||A(t)||^ < C e~ ilTH for all t > and 
any s > 0, see (|16p . so it remains to estimate r(t, x, z). We start from equation (|32|) which, in 
view of the global bound obtained in Theorem 12.11 and the estimate above for A, implies 

^l|Vr(t)||| 2 + i||Ar(t)||| 2 < C 1 \\Vrmb\\^(f)\\h+C 2 e-^ t , (51) 

for some constants Ci,C 2 > (depending on the initial data). Fix < \i < 2ir 2 and let 
f(t) = e M *||Vr(i)|| 2 2 . Using (|5ip and the Poincare inequality ||Ar|| i 2 > 2-7r||Vr|| i 2, we find 

/'(*) < e^^\\Vr(t)\\l 2 - ±\\Ar(t)\\h + C x || Vr(i) ||| 2 1| An(i) ||| 2 + C 2 e^ H ) 

< d/^IIA^Hl, + C 2 e-^ 2 -^ . (52) 

Since / °° || An(t)|| 2 2 dt < 00 by (08]), it follows from §52$ that f(t) < C 3 for all t > 0, hence 
||Vr(t)|| i2 < C 3 e~^ 2 for some C3 > 0. As H^ll^fi ~ ||V£t||£2 < ||Vr||^2 + ||VA||£ 2 , this proves 
that u(t) converges exponentially to zero in i/ 1 (B) 3 as t — > 00. The decay rate we have obtained 
so far is not optimal, but it is sufficient to conclude the proof of Theorem 11.11 

To get the optimal decay rate, the simplest solution is to go back to equation ([i~3]) satisfied 
by u. Using straightforward estimates to bound the nonlinear terms, we arrive at the differential 
inequality 



d_ 

dt 1 



|Vu(t)|| 2 2 < -2||Au(t)|| 2 2 + [ Au(t) ■ N 3 (t)dxdz 

Jo 

< -2||An(t)||i 2 + C||A^)|| L2 ||V^)|| i2 (||VS(t)|| L 3 + ||V^(t)|| L 3) , (53) 
where C > is a universal constant. Now we observe that 

POO 

2 1 IIT7,t7/'-AI|2 \ jj / ri / I'll V7„~: /-A l|2 i II V7„-; _ f+MI 2 i II l|2 



(\\Vu(t)\\i s + \\Vu(t)\\h)M < C (||Vfi(t)||is + ||Ve 3 (t)|li» + ||w3(t)||is)di < 00. (54) 

J 

For u andn 3 , this claim follows pgjl . (09]) . because ||Vu|| 2 3 < C||Au|| 2 2 < C(|| Ar|| 2 2 + 1| AA|| 2 2 ), 
and HV^all 2 ^ < CIlVusll^HAuall 2 ^ 3 < C(ll V^H 2 ^ + ||An3||| 2 ). On the other hand, the decay 
rates established in Section 13.31 below will show that ||a> 3 (i)|| L 3 = 0(t" 2/3 ) as t — > 00, so that 
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(p3|) holds. Combining (fo^j) . (f5Tj) . and using the Poincare inequality ||An|| L 2 > 27r||Vn|| L 2, we 
easily obtain 

sup e M *||V«(t)|| i 2 < oo , for any fi < 4ir 2 . (55) 

Note, however, that the linear decay rate \x = An 2 cannot be reached by this argument, because 
J °° ||Vu(t)||£3 dt = +00 in general. 

For later use, we mention that similar decay estimates can also be obtained for ||An|| i 2, by 
differentiating (|13p and repeating the same arguments. We thus obtain 

sup e At *||An(t)|| L 2 < 00 , for any fi < 4ir 2 . (56) 
t>i 



3.2 Evanescence of u 3 

We next consider the third component of the two-dimensional velocity u, which according to 
([HI) satisfies the evolution equation 

d t u 3 + (u h ■ V)u 3 + N x = Au 3 , (57) 

where Uh = (ui,U2) t - The inhomogeneous term Ni in (I57p is clearly negligible for large times, 
because ||Aq|| L 2 < || |u| | Vu\ \\ L 2 < C\\ 

Au||^ 2 so that / °° ||iVi(t)|| L 2 dt < 00. By Duhamel's 
formula, the solution of (|57|) can be represented as 

u 3 {t) = Su{t,t )u 3 {t ) - [ S u (t, a)JVi(a) ds , t > t > , (58) 

where Su(t,to) is the two-parameter evolution operator associated to the linear convection- 
diffusion equation dtf + (uh • V)/ = A/ in R 2 . As is well-known [16j [2], the operator Su can be 
expressed by an integral formula 

(Su(t,t )f)(x) = / r u (t,x;to,xo)f(xo)dx , t > t > , 

where the kernel Tu(t, x; to, xo) has the following properties: 

i) For any (3 € (0, 1) there exists Cp > such that 

C I \ 2 

< T n (t,x;t ,x ) < j-L- exp(-/? ^ ~_ X °' ) , (59) 

for all t > t > and all x, x G R 2 . 

ii) For any t > to > and any x, 10 £ K 2 , one has 

/ T u (t,x;t ,x )dx = 1 , / T a (t,x;to,xo)dxo = 1 . (60) 

It is very important to note that estimate ([59]) holds uniformly for all t > to, with a constant Cp 
which is independent of time. This is because lj 3 = d\U2 — d^u\ is uniformly bounded in L 1 (M 2 ), 
see PS]. It follows in particular from (f59|) . (f60|) that ||5u(t, to)/IU 2 — II/IIl 2 f° r an * > *0i an d 
that Su(t,to) satisfies similar LP-L q estimates as the heat semigroup e(* -t °) A . 

We claim that the solution u 3 (t) of (|57p converges to zero in L 2 (R 2 ) as t -> 00. To prove 
that, fix any e > 0, and take to > sufficiently large so that f t °° ||JVi(s)||ia ds < e. Then 

/■£ poo 

/ 5 (t,s)JVi(s)da < / ||iVi(s)|| L 2 ds < e , for all t > t , 
Jto L2 Jt 



15 



hence in the right-hand side of ((55]) it is sufficient to bound the first term v(t) = Su(t,to)u 3 (to). 
Since 1*3 (to) G L 2 (R 2 ), we can decompose 2/3 (to) = «i + v 2 with t>i G L X (R 2 ) n L 2 (M 2 ) and 
11^2 Ilia < £• Then v(t) = t>i(t) + «2(t) with 

C 

IM*)lli 2 = II^MoMHi 2 < Q._ tQ y /2 1HU 2 7^ 0. 
and ||«2(t)||ia = ||<5u(t, to)^2 Ilia < 11^2 Hi 2 — £ - Thus, if t > to is sufficiently large, we have 
||«s(t)||ra < Ih(*)l|i 2 + lh(t)||i2+ / 11^(5)11^2 ds < 3e , 

J to 

which proves the claim. 

On the other hand, we know from (|48p that ||Vu3(t)||^ 2 dt < 00, hence there exists a 
sequence t n —>■ 00 such that ||V« 3 (t n )|| 2 2 — ► as n —>■ 00. In view of ((29]) . we have for each n: 



sup||Vu 3 (t)||i 2 < ||Vw 3 (tn)|l22+C / (||Vu3( S )||t 2 ||a>3(i)||i2 + ||V«( S )||i 2 ||Au( S )|| i2 )d S , 

and the right-hand side converges to zero as n — > 00. This shows that ||Vu 3 (t)||£ 2 — > as t — > 00, 
and we have therefore proved that u 3 (t) converges to zero in /^(M 2 ) as t — > 00. 

3.3 Diffusive estimates for u) 3 

We now turn our attention to the third component of the two-dimensional vorticity Q, which 
evolves according to (JEJ): 

d t LU 3 + (u h • V)lu 3 + N 2 = Au 3 . (61) 

By (gHJ, there exists C 4 > such that ||u> 3 (i)|| L i + ||o; 3 (t)||i2 < C 4 for all t > 0. To obtain 
sharper estimates, including decay rates in time, we use a standard method that goes back to 
Nash, see [5j. By the Gagliardo-Nirenberg inequality, there exists C > such that ||u) 3 ||? 2 < 
C 1 1 cX>3 || XZcZ> 3 1 1 , hence ||w3|| 2 2 < CC^Vo^H^. Inserting this bound into ((30]) . we obtain 

^\\Co 3 (t)\\h < -C5||^3(t)||| 2 +8|||n(t)||Vn(t)||| L2 , (62) 

where C5 = (CC^) -2 . Since || |tt(t)| |Vu(t)| ||^2 decays exponentially to zero as t — > 00, it follows 
from ((Ml) that 

sup(l + t)||w 3 (t)|| 2 2 = C 6 < 00 . (63) 



A similar argument can be used to estimate ||Vw 3 ||£ 2 . Prom (|6ip we have 

|Vlj 3 || 2 2 = - / |Au> 3 | 2 dx + / (Au! 3 )(u h -V)0 3 dx+ / (A^A^dx. 



Id 

2dt' 



Integrating by parts and using the fact that HVu^H^ = ||a^3 1| ^2 , we find 

(Au} 3 )(u h • V)w 3 dx < || I VcD 3 1 |V%| |Va> 3 ||| L i < ||Va>3||| 4 W^Wlh 

< C 1 1 A.cD 3 1 1 ^2 1| Vo> 3 ll^a ||cD 3 ll^a < C|| A.d> 3 ||^4 2 1|^3||^2 2 



/ 



hence 



^ t \\VQ 3 (t)\\h < -||A^3(t)||i2+C(||a)3(t)|||2 + ||iV 2 (t)||| 2 ). 
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As \\V0 3 \\ 2 L2 < || ^3||l2 ||Aa>3|| L 2 ^^(l + ^-^llA^II L 2 , we conclude that 

^ t \\VCo 3 (t)\\h < -C 6 - 1 (l + t)||V^(t)||l 2 +C(||^(t)|| 6 L2 + ||iV 2 (t)||| 2 ) . (64) 

Now, since ||w3(i)||| 2 < ^(l + t) -3 , and since ||iV2(£)||^2 decays exponentially to zero as t — > oo, 
the differential inequality (fMj) implies that ||Vu)3(i)|| 2 2 decreases at least like t~ 2 as t — > oo. 
Taking into account the fact that 0)3(0) G L 2 (R 2 ), we arrive at 

supt(l + t)||Vw 3 (t)||^i = C 7 < 00 . (65) 
3.4 Compactness of the rescaled solution 

To show that the solution Lds(t, x) of (foTj) converges to Oseen's vortex as t — > 00, it is convenient 
to introduce self-similar variables. Following [3 [8], we define 

fls(*, x ) = YTt w ( log(1 + 1] ' 7T+j) ' (66) 

U h (t,x) = J—, r I ]og( I + / ) . 

We also denote 



Then the rescaled vorticity u>(r, £) satisfies the equation 

d T w + • V^w + N 2 = Aj:w + -(£ • V 5 )u; + iu , (67) 

where iV^r, £) = e 2T N2(e T — 1, £ e T / 2 ), and v(t, £) coincides with the two-dimensional velocity 
field obtained from u>(r, £) via the Biot-Savart law (|8ip . It is clear that 

OO /'OO /*oo 

||iV 2 (r)|| L idr = / e r ||iV 2 (e T -l)|| L idT = / \\N 2 (t)\\ L i dt < 00 , 
io io 

hence the term iV 2 (T, £) in (|67[) will be negligible for large times. The solution of (I67p can be 
represented as 

u>(r) = S v (t,t )w(t ) - / Su(r,s)JV 2 (s)ds , r > t > , (68) 

where in analogy with (|58p we denote by S^t, to) the two-parameter evolution operator asso- 
ciated to the linear equation d T w + (v ■ V)w = Aw + g(£ • V)w + w (note that S v depends on 
the velocity field v(t,£), which is considered here as given). Using the same notations as in 
Section 13.21 we find that 

(S v (T,T )f)(O = I e T T u (e T -l,^e T / 2 ;e T0 -l^oe T °/ 2 )f(io)dCo ■ (69) 
Jr 2 

The aim of this paragraph is to prove the following basic result: 
Lemma 3.1 The solution {w(t)} t >q of (I67p is relatively compact in L 1 (M 2 ). 
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Proof. By construction w G C°([0, oo), L^R 2 )) and |Mt)|| l i < C 4 for all r > 0. To prove 
compactness, we use the Riesz criterion [18] and proceed in two steps: 
i) We first show that 

sup/ Kr, £)|d£- ► 0. (70) 

T>0J\£\>R -R^oo 

Indeed, fix e > and take tq > large enough so that \\N2(t)\\ l i dr < e/2. Then choose 
i?i > large enough so that 



sup / \w(t, £)| d£ < e . 
Te[o,T ]J\Z\>Ri 



This is clearly possible, because the finite-time trajectory {w(t) | < t < tq} is compact in 
L 1 (R 2 ). For r > t the solution of §3) can be represented as in (|68p . where the second term in 
the right-hand side satisfies 



TO 



S v (T,s)N 2 (s)ds < / \\N 2 (s)\\ L ids < e/2 . 

L J TO 



As for the first term w\{t) = S v (t,tq)w(tq), it can be estimated by a direct calculation, using 
the representation formula ([69]) and the bounds (f59j) on the kernel IV Proceeding exactly as in 
the proof of [HJ Lemma 2.5], one finds R 2 > such that 



sup / \wi(T,g)\d£ < 

T>T J\f\>R2 



i«i(r.ai<i; < " 

T>T J\£\>R 2 

If we now choose R = max(i?i, R 2 ), we see that Jiti>^ \w(t, £)| d£ < e for all r > 0, which proves 



(\>R 

UD. 



ii) Our second task is to verify that 

sup sup / \w(T,£-rj) - w(r,£)|d£ — — > . (71) 

r>0 \ V \<SJM. 2 

By compactness of the finite-time trajectory, it is sufficient to check (|71h for r > 1. Using the 
definitions (l66l) and the bound (I65T) established in Section 13.31 we find 



sup || Vio(r) || = C*8 < oo 

T>1 



Fix e > 0. By the first step, there exists R > 1 such that 

sup / \w(r,g)\ d£ < § . 

r>l7|$|>i?-l 3 

Take 5 G (0, 1] such that C^tt 1 / 2 ^ + 1) < e/3. If r? G M 2 satisfies |t/| < 5, we have 

/ \w(T,Z-ri)-w(T,Z)\6Z < 2 / Kr,0|de<^. 
J|g|>_R ^ICI>R-i d 

On the other hand, by Fubini's theorem and Holder's inequality, 

/ \w(t, £ — rj) — w(t, £)\ d£ < / \r] ■ Vu>(r, £ — rrf)\ dr d£ 

J\£\<R J\£\<rJo 



< \ v \ / \vw(T,o\az < CsH^iR + i) < £ - 
J\i\<R+i 6 



hence J R2 \w(t,£ — rj) — u>(t, £)| d£ < e for all r > 1 whenever |r/| < 5. This proves (|71|). By the 
Riesz criterion, (|70|) and (fTT1) together imply that the trajectory {w(t)} t >q is relatively compact 
in LW). ■ 
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3.5 Determination of the a>limit set 

We know from Lemma f3.1l that the solution {w(t)} t >q of (j67|) lies in a compact subset of L 1 (M 2 ) 
Let f^oo be the cj-limit set of this solution, namely 



n. 



\ Woo G L X (R 2 ) 3r n — ► oo such that w(r n ) ► Woo > . 

I n— too ) 



Since L 2 w(t, £) d£ = J R2 W3(e T — 1, x) dx = a for all r > 0, where a is given by (JH]) , it is clear 
that 

(72) 



/ Woo(0 d£ = a , for all G Ooo . 
Jr 2 



Our goal is to show that fioo = {ag}, where g(£) = 1 e l£l 2 / 4 . This will imply that 

||io(r) — a^||xi — > as r — > oo, which is equivalent to (fTOj) . 

Let $(r) denote the semiflow defined by the limiting equation 

<9 T u> + v ■ V^w = A^w -\ — £ • V^w + w , (73) 



where v is the velocity field obtained from w via the Biot-Savart law (I8ip . Note that (|73D is 
just the ordinary two-dimensional vorticity equation expressed in self-similar variables. We shall 
prove that the w-limit set of the solution w(t) of (|67|) is totally invariant under the evolution 
defined by (l73]k 



Lemma 3.2 The uj -limit set Qoo satisfies ^(t)^^ = JIqo /or all r > 0. 

Using [8, Proposition 3.5], we deduce that C {a'g \ a' G R}, hence = {ay} in view 
of (I72p . This is the desired result, which completes the proof of Theorem ll.il 



Proof of Lemma l3.2l Let {S(t)} t >q denote the Co-semigroup generated by the Fokker-Planck 
operator A + g£ • V + 1, see [7J. If u; G L 1 (R 2 ), then for any p G [1, oo) we have the following 
estimates: 

\\S(r)w\\ LP < , \\VS(r)w\\ LP < ° M }^ , r > , (74) 

47ra(r) v cl(t) 2 » 

where a(r) = 1 — e _T . Moreover ||S(t)u;||£p < e r ^ 1- ^ for all r > if w G L P (R 2 ). 

Let t^oo G fioo, and take a sequence r n — > oo such that ||u;(r n ) — it/oo||z,i -> as n -> oo. 
Since the trajectory {w(r)} T >o is bounded in L 2 (R 2 ) by (|63|) . (|66|) . we have G L 2 (R 2 ) and 
(up to extracting a subsequence) we can assume that ||to(r n ) — u>oo||lp — > as n — > oo for any 
p G [1, 2). For each n G N, let io n (r) = u;(t + r n ) and w n (r) = u(r + r n ). Then w n {r) satisfies 
the integral equation 

w n (r) = S{r)w(T n ) -J S{t- s) (v n (s) • Vw n {s) + N 2 (T n + s)) ds . (75) 

On the other hand, if we denote w{t) = ^(r)^^, we have 

w(t) = Strjwco - / S(t — s)v(s) ■ Vu)(s) ds . (76) 



19 



Subtracting ([To]) from ([73]) and using the bounds ([71]) on the semigroup S(t), we obtain for any 
p€[l,2): 

||u>n(r) - w(t)||lp < e T(1_ p ) ||w(r n ) - WooIIlp + / — ||iV 2 + s)\\ L i ds 

•>Q a(r — s) p 

+ / szr(lkn(«)llz*/8 + l|w>(s)llL4/3)||«; n (s) -u)(s)|| L 4/ 3 ds . (77) 

JO a ( r _ s )2 p 

Here we have used the fact that S(r)v ■ Vw = S(r)V ■ (vw) = e~ T//2 V • S(t)(vw), and the bound 
1 1 1 1 Li < IMIl 4 IMIl 4 / 3 — C|l w llx,4/3 which holds in view of Proposition IA.ll We first choose 
p = 4/3 and consider equation (|77p for r in some compact interval [0,T]. The first line in the 
right-hand side converges uniformly to zero as n — > oo, and in the second line we know that 
||w n (s)|| L 4/3 + ||w(s)|| L 4/3 is uniformly bounded for all n G N and all r G [0,T]. Thus it follows 
from Gronwall's lemma [llj that 

sup \\w n (T) - w{t)\\ l a/ 3 >0. (78) 

tG[0,T] ' n ^°° 

Setting now p = 1 in (fTT|) and using ([751) . we conclude that ||u; n (r) — u)(t)||£1 — ► as n — > oo, 
for all r G [0,T]. In other words w(t + r n ) converges to ^(t)^^ as n — > oo, which means that 
^(t)woq G f^oo for all r G [0, T]. Since T > was arbitrary, we have shown that $(T)f2oo C f^oo 
for all r > 0. 

To prove the converse inclusion, we fix r > and take again Woo G S7oo. If w(r n ) — ► i^oo in 
L X (IR 2 ) as n — > oo, then after extracting a subsequence we can assume that w{r n — r) converges 
as n — > oo to some wq G 0^. Using exactly the same arguments as before, we can prove that 
Woo = &(t)wq. This shows that Qoo C $(r)Ooo, for any r > 0. ■ 



A Appendix : The Biot-Savart Law in R 2 x T 1 

In this appendix we give explicit formulas for the Biot-Savart law in the domain D = l 2 xT 1 , 
and we collect a few estimates for the velocity field u in terms of the vorticity uj which are used 
throughout the paper. All these results are well-known (see [19J ) and are reproduced here for 
the reader's convenience. 

Let u : E> -> M 3 be a divergence- free velocity field, and denote by uj = curlu the associated 
vorticity field. As is explained in the introduction, it is convenient to decompose 

u(x, z) = u(x) + u(x, z) , u(x, z) = u)(x) + uf(x, z) , x G M 2 , z G T 1 , 

where u = Qu, Q = Quj, and Q is the vertical average operator defined by (UJ). Then it is 
straightforward to verify that uj = curlu and uj = curlu. Moreover, the four vector fields u, 
u, uj, uj are all divergence-free. Thus we can consider separately the Biot-Savart law for the 
two-dimensional part (u, uj) and for the three-dimensional fluctuation (u,uj). 



A.l The Biot-Savart law for (u,u). 

Since the vector fields u, uj do not depend on the vertical variable z, the relations divu = and 
curl u = uj can be written in the following equivalent form: 

(a) { ? = 9 T ' (6) ( ~ 9 f l = ? ' (79) 

\ uj 2 = -diu 3 , w \ ami + d 2 u 2 = . v 1 
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To solve the first system (a), we observe that A113 = d^u)\ — d\0J2 and we use the fundamental 
solution of the Laplacian operator in M 2 . After integrating by parts, we obtain 



2vr 7 M 2 |x - y\ 2 \U2 



On the other hand, the solution of system (b) is just the ordinary Biot-Savart law in 



t>2. 



Here, if x = (x\,X2) G R 2 , we denote x 1 - = (— x 2 ,xi). In particular, we see from (|8ip that 
the horizontal part -u^ = («i,U2) of the velocity field -u can be reconstructed from the third 
component of the vorticity <D, an observation that is used many times in the previous sections. 

In both formulas (180)) and (181 j) . the velocity field is expressed in terms of the vorticity 
through a convolution with a singular integral kernel, which is homogeneous of degree —1. Thus 
we can apply the classical Hardy-Littlewood-Sobolev inequality [H] to both cases, and obtain 
the following result: 

Proposition A.l Let u be the velocity field obtained from uj via the Biot-Savart law A8QI) . (|81[) . 
Assume that 1 < p <2 < q < 00 and \ = \~\- If & 6 L P (M 2 ) 3 , i/ten u G L 9 (M 2 ) 3 , and i/iere 
exists a constant C > ( depending only on p) such that 

IMIl<2(r2) < C||a)|| L p 



Moreover, using C alder on- Zygmund's theory, one can show that ||Vu||lp < C||aj||_c,p for 
1 < p < 00. In the particular case p = 2, we even have [|Vu[|j^a = ||cD[|^2. 

A. 2 The Biot-Savart law for (u,u>) 

The relation between u and u> is most conveniently expressed in Fourier variables. Using the 
same notations as in (1171) . we can write 



u(x,z)=( V^^e'^M* Co(x,z)=[ J2 Co n (k) e^* +2 ™* > ^ . (82) 

Observe that the sums here are taken over n G Z* = Z \ {0}, because u and ui have zero average 
with respect to the vertical variable. Since divu = and curlii = u>, we have — Au = curio), 
hence 

. / —2irin ik2 

<( k ) = i 2 + 47r2n2 ( 2 ™n -ih I u n (k) , n G Z* , kR'. (.S3) 







Since n 7^ in (j83j) . it follows that ||tt||^+i < C||tD||fl-« for any s > 0, see (fT9l) . In particular, 
taking s = and using the Sobolev embedding i? x (B) L 9 (B) for q G [2, 6], we obtain: 

Proposition A. 2 Let u be the velocity field obtained from Q via the Biot-Savart law (I83D . If 

to G L 2 (ID), then u G L 9 (P) for any q G [2, 6], and there exists C > (depending only on q) such 
that 

\\u\\lvVB) < C||^IIl2(o) . 
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A. 3 The Leray projector 



In the Fourier variables defined by (I17p . (|18p . the Leray projector P has the following simple 
expression 

c/)-w = /-(*) + lk r-'!tL ( ■ where « = (2I) e R3 • < 84 > 

Clearly P commutes with the vertical average operator Q, which satisfies (Qf) n (k) = fn(k)6 n ,o- 
If / = Qf, we see from (|84|) that e 3 ■ (P/) = e 3 • /. In other words, the Leray projector P acts 
trivially on the third component of z-independent vector fields. 



B Appendix: Dispersive estimates 

This final section is devoted to the proof of Proposition 12.31 The arguments here follow closely 
the analysis of [H Chap. 5], and were already published in [2U] in a slightly different form. 

Proof of proposition [2~T3t If u(t, x, z) is a divergence-free solution of the linear Rossby equation 
(J9j, we first observe that the Fourier transform u n (t, k), which is defined as in (|82p . satisfies 

d t u n (t,k) + M%(k)u n (t,k) = , k el 2 , n £ Z* , 

where M^{k) is the 3x3 matrix defined by 



K\k) = (\kf + W)l + |fe|2 ^ 2n2 ( 2™ -i kl I . (85) 



—2irin ik<i 
Tin —ik 
-ik2 iki 



Indeed, the first term in (|85l) corresponds to — (Att)„(/c) = (|/c| 2 + Air 2 n 2 )u n (k). On the other 
hand, if u> = curlu, we have from (|83p 



~ /7\ ^ / 7 \ ~ /,n 6 A f ik \ 

u n {k) = ?Au n (fc) , u„(fe) = , where £ = I ^ ) . 

It follows that 

e 3 A u n (k) = — p e 3 A (£ A £ n (fc)) = ((e 3 • w n (fc))£ - (e 3 • £)u) n (/V)) . 

The last member is the sum of two terms, one of which is proportional to £ (gradient term) and 
the other orthogonal to £ (divergence-free term). Thus 

-P(e 3 Au n (/c)) = -p (e 3 • t)u n (k) = -r^T t A > 

which gives the second term in (|85p. 

As is easily verified, the eigenvalues of M^(k) are |£| 2 and |£| 2 ± zfir/, where 

2-7rn 



1^1 = |£(*,n)| = vlfHiV , and 7/ = V (k,n) = — . (86) 

^J\k\ 2 + 47r 2 n 

Moreover, the eigenvector corresponding to |£| 2 is proportional to £, whereas the normalized 
eigenvectors w^(k) corresponding to |£| 2 ± ifirj are orthogonal to £. Since u is divergence-free, 
we can forget about the first eigenvector, and we obtain the representation formula 

u n (t,k) = e-^^^iul^^tik)) + e u ^(u° n (k),w~(k))) , t>0, (87) 
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where u^(k) = u n (0, k) and (-, •) denotes the usual scalar product in C 3 . 

To estimate the norm of u in the space L 1 (R + , L°°(D)), we proceed as in [3] . Using standard 
approximation arguments, it is easy to show that 

II'^IIl 1 (R+,L 00 (D)) = SU P < U, (f) >z,2( R+)i 2( D )) , 

4>&£ 

where £ = {(/>£ Co°(D) | ||0||z/x>(r +i .z,i(bi)) < !}■ By the Parseval relation, we thus have 

/'OO r 

\\u\\LHR + ,L°°m) = sup / / V u n {t,k) 4> n (t, k) dkdt , 
.pes Jo Jw^ 

where </> n (i, k) denotes of course the Fourier transform of <p(t, •). The idea is now to replace 
into (I88p . and to estimate the right-hand side. Before doing that, we recall that the initial data 
Un(k) were assumed to vanish outside a finite ball Br in Fourier space, see Proposition 12.31 In 
view of (|87j) . the same property holds for u n (t,k) for all t > 0. Thus u n (t,k) = ip n (k)u n (t,k), 
where 

Mk) = (1-Wx(^^?) , fceM 2 , neZ. (89) 



Here x is as i n (I24p . and 5 is the Kronecker symbol. 

Given any A > and any B £ K, we denote by K[A, B] G C°°(D) the function defined in 
Fourier variables by 

K[A, B] n (k) = — e -m 2 +iBr, j fcgf 2 j n£ Z, (90) 

27T 

where |£| and r\ are as in (I86p . The following estimate will be crucial: 

Lemma B.l For any R > i/iere exists Cr > swc/i i/jai, /or any ^4 > and any B G R, £/ie 
function K[A,B] G C°°(B) de/med (JSDJ satisfies 



\K[A,B]\\ L oo m < C R 



e~^ A 
\f\B\ 



We postpone the proof of this lemma and first conclude the proof of Proposition 12.31 After 
replacing (|87|) into (|88|) , we have to estimate for each (j) G £ the quantity M_|_ + M_ , where 

Since the eigenvectors u^(fc) are normalized, the Cauchy-Schwarz inequality and the Parseval 
relation imply that \M±\ < ||tto||z, 2 (D)-^±> where 

p I roo 2 

N ±= E / e-^^M^^JtJ^dt dk 

r roo poo 

= e^ t+s ^ 2± ^- t ^Mk) 2 Mt,k)Ms,k)dtdsdk 

/*oo /"OO 

= < K[t + s,±n(s-t)]*cf>{s, •),#,•) >z»(D) dids . 

Jo Jo 
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In the last line we have used the definition (|9Up of K[A, B] and the Parseval relation again. Now, 
since <j) £ S, it follows from Young's inequality that 

| < k[a,b]*</>(s,-), ^(v) >l2(») I < \\k[a,b]\\ L!X(p) II^(*,0IIli(d) U(s,-)\\v{b) 

< \\K[A,B]\\ Laam . 
Thus, setting A = t + s, B = ±f2(s — t), we obtain from Lemma IB, II 

pOO fOO -47T 2 (t+s) /"I 

N 2 ± <C R / dtds < Cr 



o 7o vW-*l " \ Q \ 1/2 ' 

Summarizing, we have shown that \M±\ < C , R|r2| _1//4 ||'Uo||L 2 (D) f° r an <f> £ £> which in turn 
implies ||w||z,i(R +) z,°°(]D))) < (7,r|0| _1 / 4 ||uo||l 2 (id>)- This concludes the proof of Proposition 12.31 ■ 

Proof of lemma IB.lt Given A > and B £ R, we have to estimate the expression 

if[A,5](x )Z ) = -L [ Y e -A|aM)| 2 +^(fc,n) V)n(A; )2 e i(fc-x + 2 7 rn z ) dfc 
47T 2 /nj2 Z — ' 

where |£|,r/ are defined in ([86}) and ip n (k) is given by ([89]) . Here again, we follow the approach 
presented in [4, Chap. 5]. As K[A, B](x, z) is a radially symmetric function of x G R 2 , we can 
assume without loss of generality that X2 = 0. Clearly, we can also suppose that B > 0. Let L 
be the first-order differential operator defined by 

L = -1 , p 1 /, — To(l + «a(fc ) n)3 jfc2 ) ' where a(k,n) = -d k2 r)(k,n) . 
1 + Ba(k, n) z 

Then L(e lBr, ( k ' n ^) = e iBr, ( k > n ) , and integrating by parts (over the variable k S R 2 ) we find 

K[A,B]((x U 0\z) = At / V e ^(M) e i(fci^i+27rn 2 ) L t/ -A|C(fc,n)| 2 ^ n ( A; )2 N ) dyfc 
4-7H /H52 V / 

where L denotes the formal adjoint of L. A direct calculation gives 

L'(e-«MP Mk? ) = (_L_ _ l(8fea) _i^_)^l«M,l* Wtf 

We have to estimate this quantity for (k, n) G £>2,r and n / 0, because ip n (k) = if (k, n) ^ ,82/? 
or n = 0. We first observe that 

\£(k,n)\ > 2-ir , and \a(k,n)\ = io^^^Lq/o ^ 



{\k\ 2 + 4ir 2 n 2 ) 3 / 2 ~ 4R 3 ' 

Moreover, there exists Cr > such that |a(/c,n)| + \dk 2 a(k,n)\ < Cr. As a consequence, we 
have 

1 \1-Ba 2 \ \a\ Cr 



l + 5a 2 ' (l + 5a 2 ) 2 ' 1 + 5a 2 ~ 1 + 5/fc 2 ' 
so that |L*( e - A l«( fc ' n )l 2 ^J n (k) 2 )\ < Cj? e _47r2A Vn(^)(l + Bk%)~^- We conclude that 

H^IA^IIloo^) < JL^ ^ |L*(e- A ^)l 2 Vn(A;) 2 )|dA ; 



-4tt 2 A 



< C R e-^ A I < C R € -^ 



which is the desired estimate. 



24 



References 

[1] A. babin, A. mahalov & B. nicolaenko, Global regularity of 3D rotating Navier-Stokes 
equations for resonant domains. Indiana Univ. Math. J. 48 (1999), 1133-1176. 

[2] E.A CARLEN & M. LOSS, Optimal smoothing and decay estimates for viscously damped 
conservation laws, with applications to the 2-D Navier-Stokes equation. A celebration of 
John F. Nash, Duke Math. J. 81 (1995), 135-157 

[3] J.Y. CHEMIN, B. desjardins, I. Gallagher & E. grenier, Anisotropy and dispersion in 
rotating fluids. Nonlinear partial differential equations and their applications. College de 
France Seminar, Vol. XIV (Paris, 1997/1998), Stud. Math. Appl., North-Holland, Amster- 
dam 31 (2002), 171-192 

[4] j.y. CHEMIN, B. DESJARDINS, I. Gallagher & E. grenier, Mathematical Analysis of 
Rotating Fluids, Oxford University Press (2006) 

[5] E. fabes & D. stroock, A new proof of Moser's parabolic Harnack inequality using the 
old ideas of Nash. Arch. Rational Mech. Anal. 96 (1986), 327-338. 

[6] H. fujita & T. kato, On the Navier-Stokes Initial value problem, Arch. Rational Mech. 
Anal. 16 (1964), 269-315 

[7] T. gall AY & C.E wayne, Invariant Manifolds and the Long-Time Asymptotics of the 
Navier-Stokes and Vorticity Equations on M 2 , Arch. Rational Mech. Anal. 163 (2002) 209- 
258 

[8] T. gall AY & C.E wayne, Global stability of vortex solutions of the two-dimensional Navier- 
Stokes equation, Comm. Math. Phys. 255 (2005), 97-129. 

[9] H.P Greenspan, The Theory of Rotating Fluids, Cambridge University Press (1968) 

[10] E. grenier & N. MASMOUDI, Ekman layers of rotating fluids, the case of well prepared 
initial data. Comm. Partial Differential Equations 22 (1997), 953-975. 

[11] D. henry, Geometric theory of semilinear parabolic equations. Lecture Notes in Mathe- 
matics 840, Springer (1981) 

[12] T. kato, Strong L p solutions of the Navier-Stokes Equation in R m , with Applications to 
Weak Solutions, Math. Z. 187 (1984), 471-480 

[13] J. leray, Essai sur le mouvement d'un liquide visqueux emplissant l'espace, Acta Matem- 
atica 63 (1933), 193-248 

[14] E. lieb & M. LOSS, Analysis. Graduate Studies in Mathematics 14, AMS, Providence, 
1997. 

[15] N. MASMOUDI, Ekman layers of rotating fluids: the case of general initial data. Comm. Pure 
Appl. Math. 53 (2000), 432-483. 

[16] H. OSADA, Diffusion processes with generators of generalized divergence form. J. Math. 
Kyoto Univ. 27 (1987), 597-619. 

[17] J. pedlosky, Geophysical Fluid Dynamics, 2nd edition, Springer Verlag, New York (1987) 



25 



[18] M. reed & B. SIMON, Methods of modern mathematical physics. IV. Analysis of operators. 
Academic Press, New York, 1978. 

[19] v. roussier-michon, Long-Time Asymptotics of Navier-Stokes and Vorticity equations in 
a three-dimensional Layer, Communications in PDEs, 29 (2004), no 9-10, 1555-1605 

[20] V. roussier-michon, Sur la stabilite des ondes spheriques et le mouvement d'un fiuide 
entre deux plaques infinies. These de l'universite de Paris-Sud, 2003. 



26 



